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In paper I in this series, we found exact expressions for the equatorial homoclinic orbits: the 
separatrix between bound and plunging, whirling and not whirling. As a companion to that physical 
space study, in this paper we paint a phase space portrait of the homoclinic orbits that includes 
exact expressions for the actions and fundamental frequencies. Additionally, we develop a reduced 
Hamiltonian description of Kerr motion that allows us to track groups of trajectories with a single 
global clock. This facilitates a variational analysis, whose stability exponents and eigenvectors 
could potentially be useful for future studies of families of black hole orbits and their associated 
gravitational waveforms. 



I. INTRODUCTION 



The transition from inspiral to plunge is a crucial 
landmark in the radiative evolution of a compact object 
falling into a supermassive black hole. A natural phys- 
ical divide, the transition is also a natural conceptual 
divide. The inspiral can be modeled as adiabatic evolu- 
tion through a sequence of Kerr geodesies [l|, d, H, 0, 0, 0| 
while the plunge is currentl y b est modeled by numerical 
relativity 0, i, i, 0, [HI, El, El- Inspiral gives way 
to plunge through an important family of scparatriccs. 
In paper I in this series [iJI, we detailed the nature of 
the separatrix between bound and plunging orbits as a 
homoclinic orbit - an orbit in the black hole spacetime 
that whirls an infinite number of times as it asymptotes 
to an unstable circle. We found exact solutions for the 
family of homoclinic trajectories and depicted them as 
the infinite limit of a sequence of zoom-whirls El. As 
a companion to that physical space picture, we analyze 
the complementary phase space picture here. 

As discussed at some length in paper I, formally, the 
homoclinic orbit lies on the intersection of the stable and 
unstable manifolds of a hyperbolic invariant set. In the 
black hole spacetime, the hyperbolic invariant set is rec- 
ognized by the more familiar tag "unstable circular or- 
bit" . To make this connection precise from the phase 
space perspective, we examine the variational equations 
- the equations governing the evolution of small dis- 
placements from the circular orbits. It is straightfor- 
ward to show that the energetically bound, unstable cir- 
cular orbits are hyperoblic; that is, they have an unstable 
eigendirection and a stable eigendirection. We then show 
that the stable and unstable eigendirections arc tangent 
to the homoclinic orbit in the local neighborhood of the 
unstable circular orbit. In other words, two of the eigen- 
solutions of the variational equations around bound un- 



stable circular orbits are local representations of the ho- 
moclinic orbit. These eigensolutions capture the qualita- 
tive and quantitative features of the separatrix discussed 
in paper I, including the azimuthal motion E3l- 

We begin by devising a reduced Hamiltonian formu- 
lation of equatorial Kerr motion that natrually admits 
comparisons of groups of trajectories against a single 
global clock. The variation of Hamilton's equations yields 
stability exponents for circular orbits that could have 
general utility, for instance, as an estimate of inspiral 
or merger timescales [isl Efl| , or in a coarse graining of 
the template space around periodic orbits E3 ■ For com- 
pleteness, we also find explicit expressions for the actions 
and the frequencies 



II. KERR HOMOCLINIC ORBITS IN PHASE 
SPACE 



Carter famously reduced the full geodesic equations of 
motion to four first order equations in space and time 
coordinates E3l ■ Despite the appeal of this accomplish- 
ment, a phase space analysis requires variation of the full 
equations of motion for both the coordinates and their 
conjugate momenta. For this reason we will not work in 
the first-order integrated system of equations, although 
we will borrow his familiar expressions. Instead, we write 
down a Hamiltonian formulation of Kerr geodesic motion 
and explicitly derive the equations of motion. 



A. Kerr Equations of Motion 
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Although written out in many places, including paper 
I E3l , to remain self-contained we include the Kerr met- 
ric in Boyer-Lindquist coordinates and geometrized units 
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(G = c=l): 
ds^ = - 



2Mr\ , , 4:Marsm'^0 , , 
at dtdif 



sin^ d ^r^ + 

^dr^ + Ed6|2 
A 



2A/a2rsiii^i 



, (1) 



where M, a denote the central black hole mass and spin 
angular momentum per unit mass, respectively, and 



S = + cos^ 6 
A = - 2AIr + (? 



(2) 



The constants of motion along Kerr geodesies are the rest 
mass of the test object, energy _E, axial angular momen- 
tum L^, and the Carter constant Q [Tsj . 

In dimensionless units, the first-order geodesic equa- 
tions are 



T.t = - — - — '— - a^E sin^ ( 



(3a) 
(3b) 

(3c) 
(3d) 



where an overdot denotes differentiation with respect to 
the particle's (dimensionless) proper time t and 



6(61) = Q -cos^e la^{l 



R{r) ^ -(1 - E'^Y + - [0^(1 - E^) + i^] 



2{aE- 



QA 



(4) 



(5) 



The four equations ([3]), though no doubt valuable in 
many contexts, do not lend themselves to a variational 
analysis. The formalism we will imploy is Hamiltonian 
and a phase space study requires not just the coordi- 
nates but also their conjugate momenta. Although we 
start from scratch with a Hamiltonian formulation of the 
dynamical equations, we will make use of the Eqs. ([S])-© 
along the way. 

As in paper I, we will restrict attention to equatorial 
orbits and defer non-equatorial motion to a future work. 
Equatorial Kerr orbits have 9 = 7r/2, ^ = 0, and Q = 0. 



B. Hamiltonian formulation 

The Hamiltonian for a relativistic non-spinning free 
particle of mass fi is [l^ 



(6) 



where the inverse metric components g"^ are functions 
of the spacetime coordinates and each pa is both a com- 
ponent of the 4-momentum one-form and the canonical 
momentum conjugate to coordinate q". 

We want to build the Hamiltonian explicitly from Eq. 
(P), and we could do so just by inserting the inverse 
metric and turning the crank. However, we can yield 
an equivalent but algebraically nicer expression for the 
Hamiltonian with far less effort. To begin, consider the 
terms in the Hamiltonian explicitly containing pr or pg: 



(7) 



Since the r, 6 portion of the metric g^,^ is diagonal, that 
block of the inverse metric is also diagonal, with g^^ = 
l/grr a-nd g^^ = 1/gee- The Pr,Pe terms in H are thus 



(8) 



The remaining terms in the Hamiltonian will be 
quadratic in the remaining momenta pt and p^ with co- 
efficients that are functions only of r and 9 (since the 
metric, and thus the inverse metric, are cyclic in the t 
and (p coordinates). The Hamiltonian can therefore be 
written as 

^^(q,P) = I (^^^pl + l(^^^pl + ^F{r,9,pt,p^) , 

(9) 

where F{r,9,pt,p^) = F{r,9, E, L) is some expression 
equivalent to g"pj -I- 2g^fptp^ + g'^'^p^,- 

Notice that the f and 9 equations of ([3|) can be recast 



as 



A 2 i? „ 

— v = U 

2E^'' 2AS 

1 2 e 

2S^'^ " 2S " " 



(10) 



Adding these equations and subtracting 1/2 from both 
sides tells us that 



A 2^ 1 2 



21]"'^ 21]"" 2AI] 2S 2 



R e 1 _ 1 



(11) 



Since H = —1/2, the left hand side must be identical to 
H . Matching to Eq. ([9]), we glean that 



AE 



so that we finally get 



A 2 1 2 + Ae 1 
" 2E^'' + 2E^'' " 2AE " 2 



(12) 



(13) 



where R and 8 are the functions in ([5]). Note that in 
dimensionless coordinates, the Hamiltonian has the same 
constant value —1/2 along any trajectory. We also used 
this form of the Hamiltonian in Appendix A of Ref. [l3l • 
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Because all dependences on E = ~pt and = Pip applied to the Hamiltonian (|13p yield equations of motion 
are locked inside R and Q and H is cyclic in t and ip, 
Hamilton's equations 

■M . dH , ^, 



i 



A 
1 



-P0 



1 



d 



2AS aL 
d_ 

2ASaB 



(i? + Ae) , 
(i? + Ae) , 



Pe = - 
P^ = 
Pt = 



A 

2S 



21] 
2E 



i? + Ae 
2AS 

i? + Ae 

2AS 



(15a) 

(15b) 

(15c) 
(15d) 



where the superscripts ' and 9 denote differentiation with 
respect to r and 9, respectively. Notice, all of the Eqs. 
(fT5|) are dynamically equivalent to Eqs. ([3]). These equa- 
tions define an 8D phase space, one axis for each of the 
4 coordinates t, r, 9, ip and their corresponding conjugate 
momenta, with r parametrizing trajectories in the space. 
The Hamiltonian (fT3|) derived above governs the evolu- 
tion of the system in this 8-dimensional phase space. 

A manifestly covariant form of Hamilton's equations, 
equivalent to p4)) . has been used in other references to 
deduce important information about individual trajecto- 
ries [ll,[il,|23|. We, however, want to describe how mul- 
tiple trajectories evolve relative to one another to locate 
stable and unstable flows in phase space, and that task 
requires tracking evolution with respect to some global 
clock. In the covariant Hamiltonian picture, the time pa- 
rameter T in (|14p flows differently on different trajectories 
and is thus not a physically viable global clock. ^ 

Coordinate time t would be a good global clock, but 
it becomes awkward to maintain the clock as a coor- 
dinate in the 8D phase space. Furthermore, all orbits 
move monotonically away from the origin along the t di- 
rection.^ Consequently, no region of finite phase volume 
contains any orbit in its entirety, and there are no recur- 
rent invariant sets.'^ The 8D space, then, is not a natural 



backdrop for the discussion of homoclinic orbits. 

Indeed, this lack of boundedness is the hallmark of 
relativistic systems, in which time itself is a coordinate. 
Luckily, we can work in a 6D space - the phase space of 
spatial coordinates and their conjugate momenta - pa- 
rameterized by coordinate time t. To do this properly, 
we work with a new Hamiltonian function, the energy 
E, that generates the flow parameterized by coordinate 
time,"^ 

dq' _ dE dpi _ dE 

'dt ~djH' m ^ ' 

For details of the phase space reduction formalism see 
Refs. [2l|. It must be stressed that we treat every E in 
the Hamiltonian (|13p as an implicit function E{q,p) of 
the spatial and pi and solve 

H{q,p,E{q,p)) = -]^ (17) 

for 

In other words, the spatial part of relativistic free 
particle motion maps to an equivalent classical problem 
for which coordinate time t is the time parameter and 
whose dynamical evolution is governed by the Hamilto- 
nian E{q,p). Such a space-time splitting, which we also 



^ Mathematically, of course, r is a perfectly fine global clock. After 
all, the Hamiltonian formalism knows nothing about relativity 
and is perfectly happy to answer physically unsensible questions 
like how equal r separations evolve with respect to "global proper 
time" . 

^ Strictly speaking, the motion is also monotonic in the i/p direction, 
but topologically identifying ip = and = 2tt compactifies 
phase space in the </? direction and thus bounds the tp motion. [Tgll 

^ Of course, every inidividual trajectory is still a trivial sort of 
invariant set. Since even in this space, the phase trajectories 



describing the orbits in paper I asymptote at r — > ±oo to those 
representing unstable circular orbits, we can still talk about their 
being homoclinic to an invariant set. Still, the language is inele- 
gant, and having to track the additional t evolution is an unwel- 
come complication. 

Simply restricting attention to the spatial 6D subspacc of the full 
8D space is not formally equivalent to using the non-covariant 
Hamiltonian. We elaborate on this in future work. 
^ Since we consider only positive energies, we keep the larger root 
in the resulting quadratic equation for E. 
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used in [17[ and a fuller discussion of which we are devel- 
oping in a coming work, is dynamically exact and involves 
no approximation. The only cost is that the accumula- 
tion of proper time r along any trajectory (for which we 
will have no need in this paper anyway) must now be 
tracked on the side as a separate function.^ 



To get the 6D equations of motion for the Kerr system, 
we could calculate E{q,p) explicitly from pT]) and then 
apply (fT6|) . Alternately, we can realize that we have to 
get the same result if we divide all the spatial equations 
in ifTS]) by i (|15dp and immediately write down 



dr 



1 A 



dt~ i"" j:^' ' 



dpr 

~dt 



de 1 1 

-dt^l'^^p'^ 

dip _1 ( Id 
lit ^1^ \2K^dL 



dpe _ 1 
dt ~ i 



A 



A_ 

2E 



1 

2S 



2AY. J 
2AI] 



(18a) 

(18b) 
(18c) 



with the caveat that, when we calculate derivatives of 
Eqs. (US]), every instance of E be treated as a function 
E{q,p) rather than as either a phase space coordinate or 
a parameter. 

This 6D phase space makes variational analysis 
straightforward: because coordinate time t is both a good 
global clock and the time parameter for (|16p . the equa- 
tions dictating the evolution in t of small separations be- 
tween trajectories at equal t can be derived just by lin- 
earizing Eqs. (jlSp . We perform that linearization now. 



C. The variational equations 

We work exclusively in the 6D phase space and intro- 
duce the following notational simplification. Because the 
distinction between g's and p's as components of vectors 
and one-forms, respectively, has to do with their behavior 
in the 4D manifold of the Kerr spacetime and not with 
their function in the phase space, where they are merely 
coordinates labeling points, we will henceforth drop the 
superscript/subscript distinction. Instead, we will refer 
to both and pi as components Xi (with a subscript) of 
a single six-dimensional coordinate vector 



Pr 



(19) 



This allows us to write Hamilton's equations in the com- 
pact form 

f ^f(X) , (20) 

where the components of f can be read off Eq. (|18p . 



Now consider an arbitrary reference trajectory X(i) in 
phase space and the vector S'X.{t) of small displacements 
from points on X(t) to points at the same coordinate time 
on neighboring phase trajectories. The first order equa- 
tions of motion for (5X(i) are the linearized full equations 
of motion (PU)) around X(t). Specifically, 



d(5X(t) 
di 



or, componentwise, 

ddXiit) _ 
It ~ 

i^,,(X)^ 



X(t) 

K(X(t)) SX{t) 



(21) 



X 



dXi 



dfi 
dXi 



fixed 
E 



dE dh 
dXj dE 



(22) 



(23) 



where the last equality stems from the caveat reagarding 
equations 

Equation (|2ip is a system of first-order linear ordinary 
differential equations whose coefficients Kij(t) depend 
implicitly on time through the solutions X(t) to (|20p . 
The solution to such a system can always be expressed 
in terms of a fundamental matrix [2^ L(i;Xo) that de- 
pends on the point Xq on the reference trajectory at 
which we define the initial displacement vector JXq and 
that satisfies 



SX{t) = L{t; Xo) (5Xo 



(24) 



where L(t ~ 0; Xo) is the identity matrix. 

The goal of variational analysis is to find L, which we 
can equivalently think of as the time evolution operator 
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for small displacements. Given the equations of motion 
(|20p . we can always calculate the matrix K, but in gen- 
eral there is no corresponding analytic expression for L. 
However, K on equatorial circular orbits is the constant 
matrix*" 



/ 



K 



1 



V 



R 
2A 








A 








\ 

±^ 

7A 



1 












2 

r 





(25) 



where R" and Q^^ are the second derivatives with respect 
to their arguments of R{r) and O{0), respectively, and 7 
is a shorthand for 



T ^(^)|r— circular orbit 

Since K is constant, L has the form^ 
L{t) = 



(26) 



(27) 



and shares its eigenvectors with K. Finding the eigenso- 
lutions of (|2ip is therefore tantamount to finding eigen- 
values and eigenvectors of K. 



D. Eigensolutions of the variational equations 



The eigenvalues A of K are solutions to 
IK -All =0 



(28) 



and come in 3 pairs of equal and opposite eigenvalues 
whose magnitude we denote as 




Afl 



1 



A„ 







(29) 

(See also 0) The eigensolutions associated with the Xg 
and A^ = eigenvalues are extremely revealing in their 
own right. Presently, however, our concern is the eigen- 
solutions associated with A,-, and we defer a complete 
discussion of the eigenvectors of K to a future work. 
The Xr may be real or imaginary depending on the sign 

of 

R" 



= 12r [1 - (1 - E^)r] - 2 [a^l - E^) + i^] 



.1/2(^2 _g^^ 8^^1/2 _ 3^2) 

r3/2 _ 3ri/2 ± 2a 



(30) 



^ Although Eq. 1251 1 can be expressed solely in terms of the black 
hole spin a and the constant radial coordinate r of the circular 
orbit, we have left it in this form for readability. 

* Considerable analytic insight into L is also possible when the 
K(t) is periodic in time t, a situation that arises when the ref- 
erence trajectory X(t) is itself periodic and which we tackle for 
Kerr orbits in a future work. 



where we have used the {E, L^) found in Ref. [23| and 
used in paper I [iJl to write R" in terms of r alone. The 
plus/minus signs indicate prograde/retrograde. On the 
unstable circular orbits of interest to us (ribco < r < 
Hsco), R" is positive and A^ is real and plotted as a func- 
tion of r for various values of a in Fig. [1] 
The (unnormalized) eigenvectors 



UL 



A, Jfl, 0, 0, T- 



A -a/— ±' 



'\/«"/2 



2r^ 



'/fF72 



(31) 
(32) 



associated with ztA^ are also real. Combining and 
P7p . each eigenvalue/eigenvector pair yields a corre- 



vs r 

(UNSTABLE BOUND Prograde Circular Orbits) 

a^O - 
a=0.1 

a=0.5 - 
a=0.8 

a=0.9 - 

\ \ a=0.95 - 

5 \ a=0.998 - 



vs r 

(UNSTABLE BOUND Retrograde Circular Orbits) 




a=0 

a=0.1 

a=0.5 

a=0.8 

a=0.9 

a=0.95 

a=0.998 



FIG. 1: The dimensionless real- valued stability exponent Xr 
(measured in units of M~^) for unstable circular orbits with 
E < 1 for various spins a. Left: Prograde orbits. Right; 
Retrograde orbits. 
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spending eigcnsolution 



C(.)g-A.t^(.) 



(33) 



to the variational equation ((2T|) . where the constants c'"'°' 
refleet where we choose to set < = 0. 



E. Relation to the homoclinic orbits 

We now build the case that in the neighborhood of 
X'^"''^(t), the linearized solutions X'"-'^{t) coincide with 
exact homoclinic solutions X'^'^(t). For simplicity, we fo- 
cus first on the unstable solution in ([55]) . which corre- 
sponds to a linearized solution 



X'"'(i) =X"''"(f) + (5X<"'(f) 



(34) 



to the full equations of motion (|20p . 

Some of the similarities between the linearized and ho- 
moclinic orbit arc self-evident. The absence of 9 and pe 
components in X<">(t) indicates that the orbit remains 
equatorial, and the identical signs on the r and pr com- 
ponents reflect the fact that small displacements from 



the circular orbit along the eigendirection run away ex- 
ponentially to larger radial positions and velocities on 
an e- folding timescale A~^. The absence of a p^ com- 
ponent in (SX*."' (t) indicates that the linearized orbit has 
the same angular momentum as X'^'"^(i). 

Less self-evident is the fact that, like the homoclinic 
orbit, the linearized orbit also has the same energy E 
as the circular orbit. To see this, note that since the 
Hamiltonian E = E(X.) is a function of the phase space 
coordinates, the energy difference SE = E™^ — i?'™ can 
be expanded as a power series in the components of (5XJ."' . 
Because the derivatives of all phase variables except 
vanish on the circular orbit and 5p^p = 0, the first order 
contribution to that expansion vanishes, 



5x' + 



dH. 



6D 



dpi 



Spi 



dpi 
' dt 



Sx 



dx^ 
'dt 



5pi 



(35) 



dip 



—Spip = 



The second order variation in the energy becomes 



SE 



(2) 



dx^dx^ 

d dpj 
~dx^~dt 

d dpr 
dr dt 



Sx'^Sx^ 



dpidpj 



SpiSpj 



,d^Ihp 

dx'-dpj 



d dx"^ 

Sx^Sx-' 4- — 

opi dt 



c c ^ ^ dx^ 



Sx^Spj 
Sx^Spj 



^ 2 _|_ d dr 
dpr dt 



^ n d dr 
''^^'d-rdt 



Srdpr 



Sr^ + Kr 



5pl + 2Kr 



SrSpr 



Using Eq. ([25]) and the fact that 

1 [r^ 

on the eigcnsolution, we find that 
SE^'^ 



(36) 



5r^ 



R" 
2A2 



6r^ fR" / 1\R" ^ 




2Krr 



R" 

2A2 



A similar result holds for X''''(f), despite the addition of 
an overall minus sign in (|36[) , since through second order 
6E depends on dp^. Continuing this process to higher 
orders is beyond the algebraic patience of the authors, 
but at least through second order in the variations, the 
linearized solutions describe orbits with the same E and 
L as the unstable circular orbit. 

The component of JX^"' (t) merits more discussion. 
The ratio 6ip/Sr is fixed, so that Sip does not merely rep- 
resent an arbitrary overall translation in (p. Instead, this 
component indicates how the phasing difference between 
the linearized orbit and the circular orbit changes as the 
radial separation between the two orbits grows. Notice 
also that since 5XJ."'(i) ^ as i ^ — oo regardless of 



7 



how c'"* is chosen, the hnearized solution describes an 
orbit that is in phase with the circular orbit in the infi- 
nite past. As discussed in paper I [3], there is a unique 
choice of phase for a homoclinic orbit that will synchro- 
nize it with the circular orbit in the infinite past. Ap- 
parently, the linearized eigensolution goes so far as to 
select the phase of the homoclinic orbit it locally approx- 
imates.^ The import is that the linearization captures 
detailed information about neighboring orbits, including 
phase information. 

Analogously, the linearized solution 



(37) 



synchronizes with the circular orbit ai t = +oo. We 
can now understand the signs of the Sip components of 
both eigenvectors. In (5XJ."'(t) it has the opposite sign as 
dr because as the displaced orbit moves to larger r, its 
dif/dt drops, and it lags the circular orbit with which it 
was synchronized at t = — oo. In 5X^'''(t), in contrast, 
5ip and Sr have the same sign: since the circular orbit 
will accumulate azimuth faster than the displaced orbit 
as it spirals in, it must begin ahead of the circular orbit 
in phase if the two arc to synchronize at t = 4-oo. 

Now, as discussed in paper I [l^ . the two linearized so- 
lutions X<"'(t) and X''''(t) do not coincide with the same 
homoclinic orbit, but rather with two homoclinic orbits 
that differ by a phase. Since circular orbits that differ by 
a phase belong to the same invariant set, we continue to 
refer to these as homoclinic and not hctcroclinic trajec- 
tories. 



F. Phase portraits 

To make the coincidence between the linearized solu- 
tions and the homoclinic orbits manifest, we examine a 
phase portrait of the homoclinic orbit and the linearized 
solutions. Again, we use the radial coordinate r along 
the homoclinic orbit as our global time parameter. The 
required expression for pr in terms of r for the homoclinic 
orbit follows from Eqs. (|15ap . The result is 



Pr{r) 



(38) 



for outbound motion and the negative of the same expres- 
sion for inbound motion. Together with the exact solu- 
tions from paper I ((38)) generates the exact phase 
curves of the homoclinic orbit. Fig. [5] overlays a homo- 
clinic orbit and the corresponding linearized orbit X'"'. 
By construction, the orbits arc coincident at t = — oo. 

For illustration, we have plotted the case a = 0.8 with 
an associated unstable circular orbit at r„ = 2.500536. 



Since both orbits are equatorial (so that 9 motion can be 
suppressed) and have the same L^, a 3D orbit in r,pr,ip 
space captures all the dynamical information, and each 
panel of Fig. [5] shows the projections of the two orbits 
into a plane. The curves in Fig. [2] are the coordinate 
separations between the homoclinic and circular orbits, 
with the various projections of the separation eigenvec- 
tors overlayed. They confirm the claim made in paper 
I [31 that the global stable and unstable manifolds of 
the circular orbits are tangent at the circular orbits to 
the local stable and unstable manifolds defined by the 
eigensolutions to the variational equations. 



G. Action-angle variables 

In an action-angle formulation [3, HH, [3l of Kerr mo- 
tion, the Hamiltonian is reformulated in terms of con- 
stant momenta Ji called actions and canonically conju- 
gate angle variables ipi that increase linearly with time 
at rates LOi. Fourier expansions of orbit functionals in 
terms of the fundamental frequencies coi are the basis of 
frequency-domain radiative evolution codes, and Ref. po| 
develops a description of the inspiral dynamics entirely 
in terms of action-angle variables. For completeness, we 
include exact expressions for the frequencies and actions 
of homoclinic orbits. 



1. Fundamental frequencies 

Because the equatorial Kerr system is two dimensional 
and integrable, every bound orbit has an associated pair 
of fundamental frequencies^*^ 



2tt 

1 
T'r 

2 f ' " dr 



dip 



dt 

dip 
dr 



dt 



(39a) 



(39b) 



f dr f 

Jrp dr 

Because their radial period is infinite, w,. = for ho- 
moclinic orbits. Homoclinic orbits also whirl an infinite 
amount as they approach their periastron r„ , so both the 
numerator and denominator of (j39bp diverge. 

However, as we show in paper I [l^ . the divergences 
in both Tr and the accumulated azimuth ip can be traced 



Of course we can have a homoclinic orbit of any phase still lino 
up with the linearized solution simply by adding an overall if 
shift to SX.i"^{t). 



Even equatorial orbits have a third frequency ijjg associated with 
small oscillations about the equatorial plane. We discuss the sig- 
nificance of these frequencies for all equatorial orbits in a separate 
work. 
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FIG. 2: Projections of the eigenvector uj,"', to which the linearized separation JX'"' is proportional, overlayed with the actual 
coordinate differences X'^'^ — X'^"^'^ in the phase space. In the /S.ip plot, we have identified — vr at the bottom of the plot and tt 
at the top. The plots, intended to be schematic, are around an unstable circular orbit at ru = 2.2 for a = 0.8. 



to specific terms of the form 



2. Actions 



2—^ tanh 
A,. 

2— tanh^ 





- r 




ru 




r 




ru 



(40) 

dip/dt{ru), the con- 



Their ratio thus converges to fly 
stant coordinate velocity of the circular orbit at r^. 

The azimutlial frequency for the homoclinic orbit and 
its associated unstable circular orbit are thus the same, 



,hc 



,hc 



0, 



I a, 



(41) 



3/2 

ru 



That allows us to make a nice statement: the stable and 
unstable circular orbits determine the lower and upper 
bounds, respectively of the w^'s of all eccentric bound 
orbits with a given Lisco < < iibco- 



Each action Ji of a bound orbit is defined by 



Ji= f pi dqi 



(42) 



where the integral is taken over the projection of the 
orbit into the qi,Pi plane. Since = is constant, 
J(p = 2Tr Lz for any orbit. The radial action is the area 
enclosed by closed {r,pr) curves like that of Fig. 



Pr{r) dr 



dr 



(43) 



For arbitrary orbits, ((43|) at best reduces to elliptic inte- 
grals, but for the homoclinic orbit, Jr can be written as 
an exact function of r„ alone. The result, derived in the 
Appendix, is 



2E^ - 1 



= 2VT^ X {ra r„) + 2-^ ^ — ^ / 



1 -£;2 



tan 



Vi?(r_)tanh"^ 



r- ra - ru 



rn - r_ r„ 



a/ R{r+) tanh 



r. r„. - r,, 



rn - r. r„ 



(44) 



r 



III. CONCLUSION 



Although the results of this paper are self-contained, 
the phase space portrait is a direct complement to the 
physical space portrait of paper I [3]. Both approaches 



identify the separatrix between bound and plunging or- 
bits with a homoclinic trajectory that whirls an infinite 
number of times on asymptotic approach to a circle. 

Although the intention was to detail a profile of the 
separatrix, the technical results of this paper could have 
further utility. In partcular, the whirling stages of tra- 




FIG. 3: Schematic of a homoclinic tangle. The curve above 
represents the repeated intersection of a homoclinic orbit of 
the perturbed system with the r,Pr phase plane. The large 
dot represents the intersection of a circular orbit in the per- 
turbed system along with the eigenvectors denoting its local 
stable and unstable manifolds. 



jectories in the vicinity of the homoelinie set might be 
modeled as variations around the circular orbit using 
the eigenvectors and eigenvalues found here. In the fu- 
ture, we aim to generalize this approach to capture orbits 
around the periodic set p7| and to move out of the equa- 
torial plane [2^, [2^ . 

Another connection that should be made in a dynam- 
ical discussion of the separatirx is its role as the divide 
between chaotic and non-chaotic behavior. The geodesic 
motion of a non-spinning test particle around a Kerr 
black hole is known to be integrable [l^. There are 
as many constants of motion as there are canoncial mo- 
menta in this Hamiltonian system and the motion can 
therefore be confined to regular tori in an action-angle 
set of coordinates. 

However, the presence of a homoclinic orbit indicates 
the Kerr system is vulnerable to chaos [13, [H, [H, HO] • 
Under perturbation, the stable and unstable manifolds 
that previously coincided along the homoclinic orbit (Fig. 
[2]) can develop transverse intersections. In other words, 
the stable and unstable manifolds do not coincide but 
rather intersect, and once they intersect, they do so an 
infinite number of times creating a homoclinic tangle, 
as in Fig. [3l The homoclinic tangle is associated with 
a fractal set of periodic orbits and marks the locus of 
chaotic behavior. Chaotic behavior has in fact already 
been found in the Kerr system for spinning test particle 
motion (28j and in the case of spinning comparable mass 
black holes 

Chaos may be dissipated by gravitational radiation 
losses [H, [SJ, [3^. However, due to the poverty of 
the approximation methods in the strong-field, there is 
no definitive resolution to the question of the survival 
versus extinction of chaos in astrophysical systems. If 
chaos does survive radiative dissipation in rapidly spin- 
ning black hole pairs, the highly non-linear character of 
black hole spacetimes could be evidenced by the destruc- 
tion of the homoclinic orbit on transition to plunge. 
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APPENDIX A: DERIVATION OF ACTION OF 
HOMOCLINIC ORBITS 

The radial action of a bound non-plunging orbit is the 
area enclosed by its projection into the r, plane, 



Jr= ipr{r)dr^2 I dr^l^ , (Al) 



where and arc the periastron and apastron, respec- 
tively, and R{r) is the function ([5]). 

For a homoclinic orbit, Tp equals r„, the radius of the 
associated unstable circular orbit, and Ta is expressible 
in terms of r„ alone [l3|- Additionally, R{r) factors into 



R{r)^{l- E'^)r{r -ruf{ra-r) 



(A2) 



with E the common energy of the homoclinic and unsta- 
ble circular orbit. The orbit independent quantity A can 
always be factored into 



A = (?' — r^)(r — ?'_) 



(A3) 



where r± = 1 ± yl — a? are the outer and inner hori- 
zons, respectively, of the central black hole. Together, 
the above allows us to write the radial action (|Aip of a 
homoclinic orbit as 



2V1 - 



^^ {r - ru)^r{ra - r) 
{r — r^)[r — r_) 



dr 



'a-r r{r - r„) 
r (r — 'r^)[r — r_) 



(A4) 



The integral in (|A4p can be done analytically. Under 
the change of variable 



r u 
, r = -T. -r,, 



+ 1 



J Ifa- r 2ra 
ar \ I = au ■ 



(1 +v^r 



the factors in (IA4I) become 



r — r_ = 



ra 




1 + u2 




u^ira - 




IH 


- V? 


u^ira - 


r+) - r+ 


H 






r_) — r_ 


1- 


\- V? 



(A6) 



and (jA4[) becomes 

r 



2^l^~W {ra-r+){ra-r_) 

2rW[u^-ul] ' ^^^^ 



du 



where 



(1 -f U^)^ [^2 — ^2 j |^u2 _ y2 



' + 2 
U 



(A8) 



'a ' u 



The integral in (|A7|) decomposes by partial fractions 
into 



(AiIi + A2I2 + A3I3 + A4X4) 



(A9) 



2Vl - E^ 
where the coefficients Ai are 

Ai^ra, A2^2 (r„ - 2) 

r-_(r„-r_) • (AlO) 

^3 = , , /I4 = — - 



(A5) and the functions li are 



1 



+ tan u 



1+1*2 



1 



tan ^ u 



_ 1 
^ " 2 



(iii 



1 1 



In 



u + u_ 



u2 — u2 2 



^In 



U -I- M+ 



(Alia) 
(Allb) 
(Allc) 
(Alld) 
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The right hand side of (jA9[) is easiest to evaluate in pieces. The first two terms give 



(AiXi + A2I2) 



ra ^ o (ra + 2r„ - 4) ( - - tan ^ Uu) 
1 + uf, V 2 / 



= - {ra - r„) + (ra + 2r„ - 4) (tan ^ — 



(A12) 



-\/ru {ra - Tu) + 2 — 



2^;^ - 1 



tan 



— 1 / ' a ' u 



To go from the first to the second line in (jA12[) . we have used tan ^(m) + tan ^(1/m) = 7r/2. To get the last line, 
we have used the fact that 



fa + 2r„ — ^2 



for homoclinic orbits, which follows from equating the cubic coefficients in equations ([5]) and (|A2 
The third term in (IXOl) is 



(A13) 



and likewise 



1 {tu - '■-) V?'- {ra - ?•-) 



2 

1 1 



In 



2 V 1 - ^2 
1 



In 



V 1 - 

1 / 



tanh — 



V 1 - 



tanh ^ 



T- Ta - r„ 



ra - r_ r,. 



(AM) 



Combining ((X9|l . dSHI, ([m]) and (jXTS]) . we find that 



Rir+) , -1 
^ ' tanh ^ 



(A15) 



= 2v/l -E-^x I -V'-n (ra-r„) + 2- 



2£;2 - 1 



-£2 



tan 



1 'a ' u 



+ , ^ < v'^(r_)tanh-i 



A/i?(r+)tanh"^ 



fa - r+ r„ 



(A16) 
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